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We consider the effects of charge ordering in NaV^Os (be- 
low Tsp) on the exchange constants and on the magnon dis- 
persion. We show that the experimentally observed splitting 
of the magnon branches along the a direction is induced by 
charge ordering. We find that one can distinguish between 
the proposed 'zig-zag' and 'in-line' patterns of charge order- 
ing. Only the zig-zag ordering is consistent with the experi- 
mental results regarding (i) the unusual intensity modulation 
observed in magnetic neutron scattering, (ii) the reduction 
in the intra- ladder exchange constant below Tsp, and (iii) 
the magnon dispersion along a. We estimate the inter-ladder 
exchange constant to be 1.01 meV = 11.7 K for T > Tsp- 
PACS numbers: 64.70.Kb, 71.27.+a, 75.10.-b, 75.30.Et 



INTRODUCTION 

The physical properties of the quasi one-dimensional 
(ID) spin-liquid compound NaA^Os are currently under 
intense investigation. NaV205 is an insulator and its 
magnetic susceptibility is that of a ID Heisenberg chain 
0. Below Tsp — 34 K a spin-gap opens j[| and the crys- 
tallographic unit-cell doubles along the chain-direction (b 
axis) ||. 

The nature of the spin-chains was originally supposed 
to be an 'in-line' arrangement of V +A spin-1/2 ions along 
the crystallographic b direction, separated by chains of 
V +b spin-zero ions Recently though, the notion of 

NaV205 as a quarter filled ladder compound has been 
proposed B , based on a redetermination of the crystal 
structure |l],|| and a subsequent LDA analysis Q. In 
this interpretation the electron spins are not localized 
at V-ions, but distributed over V-0(1)-V molecules gj, 
the rungs of the constituent V-0 ladders running along b. 
This interpretation has found support by recent NMR |7]] 
and Raman Jsj measurements, though far-infrared optics 
have been interpreted as evidence for charged magnons 

B- 

The nature of the state below T$p is not yet deter- 
mined. Isobe and Ueda jjj originally proposed a standard 
spin-Peierls scenario. The detection of two inequivalent 
V-ions below Tgp by NMR |?J indicates charge order- 
ing in addition to the observed formation of a singlet 
ground-state. The pattern of the charge ordering is cur- 
rently debated, several authors have proposed 'in-line' 
and 'zig-zag' patterns, based on electronic considerations 

IMi- 

The occurence of two well defined magnon-branches 



in NaV2C>5 along the a direction (perpendicular to the 
chains), as measured by neutron scattering Jig ], has come 
to a surprise and is yet unexplained. Here we propose 
that the observed splitting of the magnon branch is di- 
rectly proportional to the square of the charge order pa- 
rameter, in the case of zig-zag ordering. In the case of 
in-line charge ordering we find only one magnon branch. 
We also observe that zig-zag/in-line charge ordering leads 
to a decrease/increase of the magnetic exchange below 
Tsp along b. Experimentally a decrease of ps 20% is ob- 
served • We finally comment on the unusual inten- 
sity modulation observed in neutron-scattering [fl5| and 
find it to be consistent/inconsistent with zig-zag/in-line 
charge ordering. 



EXCHANGE COUPLINGS 

The exchange couplings in the dimerized state (T < 
Tsp) might differ from those in the homogeneous state 
(T > Tsp), due to the charge ordering (see Fig. |l]). The 
sign and the magnitude of these changes depend on the 
pattern of the charge ordering. We consider first the 'zig- 
zag' charge ordering illustrated in Fig. [l|. We define A c to 
be the charge order parameter. (1± A c )/2 are the charge 
occupation factors of the two inequivalent V-ions. A c = 1 
would correspond to complete charge disproportionation, 

2 T/+4.5 _> y+5 + 

Following the standard derivation of the exchange con- 
stant Jjj between two V-sites i and j by perturba- 
tion theory, one finds immediately that Jjy is linearly 
dependenent on the respective site occupation factors 
Pi/j — (1 ± A c )/2, i.e. Jy ~ PiPj- The dominant mag- 
netic coupling in NaV20s is between the rungs of the 
constituent ladders along b and has contributions from 
bond-exchange processes ~ (PiPj + PkPi) and from ring 
exchange ~ \/piPjPkpli where i,j,k,l denote the four 

41. We denote with J and Joo 



y/PiPjPkPl, 

sites of the two rungs [[ 
the exchange constants between rungs along b in the low- 
and high-temperature phase respectively. The zig-zag 
charge ordering (see Fig. modifies the exchange as 
J = Joo(l-A^) g. Note that the experimental values for 
J re 441 K (T < T SP ) g|[l6j and for J x ps 529 - 560 K 
(T > T S p) 0,0 indicate A^ pa 0.2, i.e. A w 0.44 jjjj. 

Another possibility of charge ordering has been dis- 
cussed in the literature, where all charges align along 6, 
leading to alternating V +A and V +5 chains |||. For this 
'in-line' charge ordering the exchange coupling between 
rungs along b has contributions ~ (1 + A^) from the 
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bond-exchange processes and ~ (1 — A^) from the ring- 
exchange. No substantial reduction of the exchange with 
charge ordering is then expected [Q. 

In Fig. [l] we illustrate a dimerization pattern for 
NaV 2 C>5 which is consistent with the observed doubling 
of the unit-cell along the a and the b direction || , leading 
to an alternation of the exchange coupling J between the 
rungs along b: Ji/ 2 = J(l ± S). The modified exchange 
integrals between rungs along a due to the zig-zag charge 
ordering for T < T SP are (see Fig. |l|): J' = J' 00 {1 + A c ) 2 , 
J" = J'ooQ- - A?) and J'" = J^(l - A c ) 2 , where is 
the exchange along a in the high-temperature phase. 
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FIG. 1. Illustration of the dimer-phase of NaV^Os with 
'zig-zag' charge ordering. V-ions with charge (1 ± A c )/2 are 
denoted by open/filled circles. The spins are in V-0(1)-V 
bonding orbitals (rungs), denoted by the short horizontal bars 
connecting pairs of nearby V-ions (not on scale, their actual 
distance is ~ a/3). The dimers are denoted by the vertical 
bars connecting the rungs. The lattice constants are a and 
b. The primitive unit-cell is spanned by n and rh and con- 
tains one dimer (= four V-ions). The orthorhombic unit-cell 
(dashed box) contains four dimers. Note that the in-line 
charge ordering would correspond to having parallel chains 
of V-ions with charge (1 ± A c )/2 respectively running along 
6. 



MAGNON DISPERSION 

We now calculate the magnon dispersion using the 
bond-operator theory As the bond-operator ex- 

pansion is a well established technique we refer to the 
literature for details and present here only the basic 
notions. The two-electron state on a dimer (involving 
the electrons on rung 1 and rung 2) can be in a spin- 
singlet or in a spin-triplet state (s'|0(i i2 )) and ijj0(i i2 )) 
(a = x, y, z) respectively), e.g. 

S t|0 (1 , 2) ) = -^(|1T,2|)-|1|,2T)) , 



4|0(i, 3) ) = ^(|lt,2i) 



V2 



|U,2T» 



Due to spin-rotational invariance one needs to consider 
only the z-component of the inter-dimer exchange S/ • S// 
(I and I' are rung-indices) which couples s^|0( 12 )) with 
^l|0(i j2 )) and vice versa. The intra-dimer matrix elements 
are 



(o (llJ0 Mf*t|o (lj3) > = + 

(0 (llJO Mf»t|0 (lj3) ) = - 



(1) 



The inter-dimer matrix elements of SfSf, defining the 
Hamiltonian in the bond-operator representation are 
products of the intra-dimer matrix elements given in Eq. 

The dimer-covering presented in Fig. [M contains four 
dimers in the orthorhombic unit-cell. It is easier to calcu- 
late the dispersion relation using the primitive unit-cell, 
which contains only one dimer. The primitive unit-cell is 
spanned by 



a b 
2' 2 



-a 36 
T'T 



(2) 



where a = 11.316 A and b = 3.611 A are the lattice 
parameters of the room-temperature structure (|^] . The 
basis vectors n and rh span a simple square lattice of 
dimers with couplings between the dimers in real space 
J' = 2 J" — J'" along ±n and — J' along ±rh. In reciprocal 
space these couplings lead to momentum dependences 
~ cos(k-n) and ~ cos(k-m). We dehne the dimensionless 
matrix element 



J ' ft a 
cos fe T — 

2Ji I x 2 



ky 2 



(3) 



J ' fu a 
cos k r — 

2Ji V 2 



3b 



2Ji 



cos (2bk y 



with k = (k x ,k y ). We employ the "linearized Holstein- 
Primakoff" approximation (LHP) and eliminate the sin- 
glet operators via the constraint that each dimer can only 
be in a singlet or in a triplet state. This procedure leads 
to the following linearized Hamiltonian for the triplet ex- 
citations fl9[| : 



■h 



E 



[(l + ek)<aV«+f 



L k,a L - 



k . a 



+ K J 



-k,r 



(4) 



The tk,a obey Bose commutation relations, [t k a , t^, a ,] = 
Sk,k'S a , a ' and they are the Fourier transform of the triplet 
Bose operators on the dimer bonds. Eq. (Q) is solved by a 
standard Bogoliubov transformation |l^], leading to the 
(threefold degenerate) magnon-dispersion 
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E{k) = JiVl + 2e k 



(5) 



Yosihama et al. jjjl have measured E(k) by inelastic neu- 
tron scattering and found a splitting for k y — n/b. Spe- 
cializing Eq. ([5]) to k y = ir/b we obtain 

E(k x , = ^Ji(Ji - J 2 ) + Ji(J' - J') cos(^^) . 

(6) 

E(k x ,ir/b) has the periodicity 47r/a. Folding back to 
the first Brillouin zone (k x S [— 7r/a, 7r/a]) we have two 
magnon branches which are degenerate at k x = and 
split at k x = 7r/a, as seen in experiment Note 
that the dispersion along a (and the splitting) is possible 
only for a finite A c ^ 0, since 

J'-J' = 2 J" J'" -J' = JU-^Al) « -Q.8J^ , (7) 

where we have used A 2 ss 0.2. The fact that the magnon 
dispersion along a is ~ A 2 for k y = 0,ir/b (in leading or- 
der in J^) is a direct consequence of the complete mutual 
frustration of n.n. ladders in NaV^Os. 



Using Eq. (§) and taking A 2 C = 0.2 we find = 
1.01 me V = 11.7 K. The anisotropy is therefore roughly 
J'oo/Joc ~ 1/45. This small value of the inter-ladder cou- 
pling is consistent with an argument by Horsch et al. || , 
that there is a partial cancellation between intermedi- 
ate singlet and triplet states in the superexchange paths 
contributing to J^. In Fig. || we compare the resulting 
magnon dispersion with the experimental result. 
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FIG. 2. Theory (solid line) and experimental (filled 
squares) results Jl5| for the magnon dispersion in NaV20s. 

The experimentally observed gap of AE = 10 meV at 
k x = determines 5: 



AE = J y/2(l + 6)6, 6 - 0.034 , 



(8) 



where we have used J = 440 K = 37.9 meV. This value 
of 6 should be considered as a lower bound to the true 
value of 6, since the LHP overestimates somewhat the gap 
for small values of the exchange alternation <5 |fT9| . The 
splitting Ei — Ei — 3.2 me V [|15[ of the magnon branch 
at k x = it /a is given by 



FIG. 3. A second possible dimer configuration for the 
dimerized phase of NaV2 0s. The symbols are the same as in 
Fig. 1. Due to the charge ordering the dimers denoted with 
(a) and (b) are not equivalent to the dimers denoted with (c) 
and (d). The R( and Rj denote the positions of the V-ions 
(not on scale). 



DISCUSSION 

We have seen that the observed splitting of the magnon 
branches along the a direction in NaV2C>5 can only be 
explained if the total coupling between n.n. ladders, 

J 'total = J' - J> = 2 J " - J "' ~ J ' does not sum U P to 
zero. As a possible source on nonzero J' ta * a i we have 

identified the charge ordering below Tsp 0. Another 
possible candidate for a nonzero J' tota i could be, in prin- 
ciple, the structural distortions occuring for T < T$p. 
There are two types of structural distortions: (i) distor- 
tions ~ 8 leading to the dimerization along b and (ii) 
distortions leading to two inequivalent V-ions. This last 
type of distortions, which are the underlying cause for the 
charge disproportionation occuring for T < T$p, could 
also alter the exchange constant directly, contributing to 



E 2 -Ei = E(w/a, tt/6) - E(Zir/a, ir/b) 



(9) 



a nonzero J' tota i- At this point we cannot rule out this 
possibility, which however would not affect the results 
presented. Let us also remark in this context, that the 
originally proposed in-line charge ordering would re- 
sult in J' = J" = J'" and thus in a vanishing J' tota i- An 
in-line charge ordering can therefore not explain the ex- 
istence of two distinct magnon branches in NaV^C^. 

We have seen that the dimerization pattern shown in 
Fig. |l| is able to explain the measured magnon dispersion 
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in detail. We note here, that another pattern (illustrated 
in Fig. ||) is also able to explain the experimental data. 
The calculations are, however, more complicated since 
there are four dimers in the unit-cell, complicating the 
Bogoliubov transformation. We have solved the resulting 
8x8 matrix numerically and found for k y = 7r jb results 
identical to those presented in Fig. ||. For a general value 
of k y the magnon dispersion arising from the two dimcr 
patterns differ slightly. Besides the two dimcr coverings 
illustrated in Fig. (|l|) and Fig. (||) dimer coverings involv- 
ing inter-ladder dimers are in principle possible. We do 
not consider them here, since the inter-ladder exchange 
is very small. 

Next we discuss the interpretation of the unusual mod- 
ulation of the intensity of the magnon branch found in 
neutron-scattering p5[. The authors of Ref. [|l5) found 
that this intensity has a period of 3n/a in k x . The au- 
thors argued that this intensity modulation could be un- 
derstood only if there is substantial magnetic scattering 
from both legs of the structural ladders S), in accor- 
dance with the notion that NaV2C>5 is a quarter-filled 
ladder compound [[|. 

The intensity modulation l m od. is proportional to the 
square of the weighted sum of the scattering phases over 
the magnetic ions in the orthorhombic unit-cell. We 
denote the positions of the V-ions by Ri , . . . R4 and 
R'^.R^ (see Fig.|). 
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ik R, 



2=1 



Pi? 



ik-R.! 



(10) 



The R; and RJ are equivalent, i.e. exp[ik • R;] = exp[ik • 
RJ]. For the zig-zag ordering the charge occupation fac- 
tors pi = (1 ± A c )/2 obey pi + p\ = 1. The charge order 
parameter A c does therefore drop out of Eq. ( |l0| ) and the 
above expression for I mo d. reduces to the one proposed in 
Ref. jlB) . The zig-zag ordering is therefore in agreement 
with experiment. I mo d. would be substantially reduced 
for the in-line charge ordering, which has pi — p\. The in- 
line charge ordering would therefore not be in accordance 
with the neutron-scattering results presented in Ref. [n5| . 



CONCLUSIONS 

We have discussed two possible types of charge order- 
ing and found that both the unusual intensity modulation 
found in magnetic neutron scattering |l5| ] and the reduc- 
tion of the exchange constant along b below T$p indicate 
'zig-zag' charge ordering. We have examined the effects 
of charge ordering upon the magnon dispersion along a 
and found, with a bond-operator study, that the splitting 
of the magnon branches found in experiment Jl5[ is a di- 
rect consequence of the charge ordering below T$p. The 
theory presented here consequently predicts the splitting 



E2 ~ E\ to vanish as T — > T$p (from below), as the 
charge-order vanishes for T > T$p P]||@] 

We would like to note that we predict a phase shift of 
7r/2 between the two magnon branches for k y — ir/b and 
k y = 0,2n/b (see Eq. (||)). For k y = ir/b the magnons 
are degenerate at k x = and split at k x = ir/a. For k y = 
0, 2n/b the magnons split at k x — and are degenerate 
at k x = n/a. 

We would like to acknowledge discussions with P. Lem- 
mens and the support of the German Science Foundation. 
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